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statement  of  the  Problem  Studied 


Increasing  the  service  lifetime  of  military  machinery  will  improve  logistics  in  the  field 
(less  down-time  due  to  maintenance  or  replacement)  and  will  free  up  resources  for 
other  pressing  needs.  To  this  end,  the  Army  would  like  to  enhance  the  lifetime  of  gun 
barrels,  which  are  limited  currently  by  steel  erosion.  As  it  is  not  possible  to  determine 
the  erosion  mechanisms  in  real  time  experimentally  in  situ,  computer  simulations  can 
help  elucidate  those  mechanisms  and  ultimately  extend  gun  barrel  lifetime  by 
suggesting  new  ways  to  stop  the  erosion. 

In  the  past  1.5  years,  a  new  quantum  mechanics  simulation  tool  was  developed  and 
interfacial  materials  chemistry  aspects  of  the  gun  tube  erosion  problem  were 
investigated.  Specifically,  ultrasoft  spin-dependent  pseudopotential  density  functional 
theory  (DFT)  was  developed,  which  provides  an  accurate  and  efficient  ab  initio 
description  of  magnetic  transition  metals  such  as  Fe  and  Cr.  DFT  was  used  to  evaluate 
the  structure  and  stability  of  a  multilayer  protective  coating  for  steel  comprised  of  MoSi2 
and  Si02,  as  an  alternative  for  the  usual  chrome  coating.  Properties  of  the  Cr/Fe 
interface  were  also  calculated  and  the  origin  of  the  stability  of  the  chrome  coating 
identified.  Pathways  for  CO  and  FI2S  adsorption,  diffusion,  and  dissociation  on  Fe  and 
Fe  alloy  surfaces  were  also  investigated,  in  order  to  understand  the  surface  chemistry  of 
common  propellant  byproducts.  These  findings  were  reported  in  five  journal 
publications,  with  two  other  publications  in  press.  Note  that  this  grant  is  of  shortened 
duration  as  it  is  the  latter  half  of  a  grant  that  originated  when  the  P.l.  was  at  UCLA 
instead  of  Princeton.  Only  those  results  published  since  the  continuation  grant  was 
initiated  at  Princeton  are  discussed  here,  and  the  reader  is  referred  to  the  earlier  final 
report  submitted  for  the  research  done  under  the  earlier  part  of  the  grant  at  UCLA. 

Summary  of  the  Most  Important  Results 

The  most  important  findings  during  tenure  of  this  grant  are  described  below,  where  the 
references  refer  to  the  publication  list  at  the  end  of  this  section.  We  first  discuss  the  new 
DFT  pseudopotential  method  we  developed  and  how  robustly  it  describes  properties  of 
magnetic  transition  metals  such  as  iron.  Then  we  review  how  DFT  was  used  to 
elucidate  the  initial  stages  of  sulfidization  and  carburization  of  steel,  and  to  quantify  the 
stability  of  a  potential  alternative  coating  for  steel.  We  verified  that  H2S  dissociates 
readily  on  Fe,  thereby  providing  both  sulfur  to  the  surface  (the  precursor  to  FeS)  and 
hydrogen  that  can  subsequently  enter  and  embrittle  steel.  We  then  went  on  to  consider 
alloy  pretreatments  of  the  Fe  surface  as  a  strategy  for  inhibiting  deleterious  surface 
chemistry,  focusing  on  FeAl  and  FeaSi.  We  found  that  neither  could  adequately  inhibit 
H2S  dissociation  reactions,  while  Si  alloying  could  help  prevent  CO  dissociation  (the 
precursor  to  carbide  formation).  Al  alloying  of  steel  has  the  opposite  effect  desired, 
namely  to  catalyze  CO  dissociation.  Given  our  failure  to  find  a  successful  alloy 
pretreatment  strategy  that  would  simultaneously  prevent  reactions  of  multiple  propellant 
byproducts  with  steel  surfaces,  we  turned  our  attention  to  protective  coating  design. 
The  origin  of  the  Cr/Fe  interfacial  strength  was  pinpointed  -  namely  strong  spin- 
correlations  and  d-d  bonding  across  the  interface  combined  with  near  perfect  lattice 
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matching  are  responsible  for  the  excellent  adhesion  of  Cr  to  Fe.  We  then  examined 
MoSi2  and  its  oxide  scale  Si02  as  a  possible  alternative  less  toxic  protective  coating  for 
steel.  We  found  that  the  ideal  strength  of  this  multilayer  coating  is  lower  but  perhaps 
competitive  to  chrome,  suggesting  it  may  be  a  useful  coating  pursue,  particularly  in  a 
composite  form,  perhaps  coupled  to  a  thermal  barrier  coating  on  top  of  the  silica  scale, 
such  as  yttria-stabilized  zirconia. 


I.  Advances  in  Spin-Dependent  Pseudopotential  Theory 

Since  the  primary  components  of  a  gun  tube  are  steel  and  chrome,  i.e.  Fe  and  Cr,  both 
of  which  are  magnetic  (Fe  is  ferromagnetic,  Cr  is  antiferromagnetic),  it  is  critical  to 
describe  their  magnetic  properties  accurately,  since  these  may  affect  other  properties 
such  as  structure  and  energetics.  The  standard  DFT  pseudopotential  method  often  fails 
to  describe  magnetic  transition  metals  because  the  conventional  pseudopotentials  are 
spin-averaged.  Our  spin-dependent  pseudopotential  theory  is  constructed  to  give  the 
correct  pseudopotential  limits  for  spin-neutral  and  fully  spin-polarized  atoms,  and  works 
between  those  limits  by  scaling  the  fully  polarized  potential  by  the  spatially  varying  net 
spin  polarization.  It  does  so  without  added  parameters,  while  providing  a  more  flexible 
description  that  self-consistently  adapts  to  the  local  environment.  During  the  term  of 
this  grant,  we  reformulated  our  earlier  incarnations  of  the  spin-dependent 
pseudopotential  theory  to  work  within  the  computationally  advantageous  ultrasoft 
pseudopotential  formalism,  showing  that  we  retain  the  previously  demonstrated 
accuracy  of  the  norm-conserving  pseudopotential  approach,  while  dramatically  lowering 
the  cost  of  the  calculation  [1].  We  showed,  for  example,  this  ultrasoft  spin-dependent 
pseudopotential  theory  reproduces  all-electron  density  functional  theory  results  for 
many  different  phases  and  magnetic  structures  of  iron,  including  the  correct  magnetic 
ground  state  and  transition  as  a  function  of  compression  and  tension  for  austenite  (fee 
Fe),  a  material  phenomenon  where  conventional  ultrasoft  pseudopotential  theory  fails. 

II.  Applications  of  Periodic  DFT  to  the  Gun  Tube  Erosion  Problem 
A.  Interaction  of  CO  and  H2S  with  Fe  and  Fe  alloy  surfaces 

H2S  is  a  common  byproduct  of  the  propellants  used  in  military  guns  and  may  be  the 
source  of  the  FeS  observed  by  electron  microscopy  in  eroded  gun  tubes,  as  well  as  the 
source  of  hydrogen  that  may  cause  hydrogen  embrittlement  of  the  steel.  We  used  DFT 
to  characterize  FI2S  and  FIS  adsorption,  diffusion,  and  dissociation  on  the  Fe(110) 
surface.  H2S  is  predicted  to  weakly  adsorb  via  a  sigma-donor  bond  on  the  short  bridge 
(SB)  or  long  bridge  (LB)  site,  with  a  binding  energy  of  no  more  than  0.50  eV.  The  H2S 
diffusion  barrier  is  predicted  to  be  small  (~0.1  eV).  By  contrast  to  FI2S,  FIS  is  predicted  to 
be  strongly  chemisorbed  on  Fe(110),  with  the  S  atom  in  the  highest  coordination  site 
and  the  FIS  bond  oriented  perpendicular  to  the  surface,  due  to  charge  transfer  from  the 
surface  to  S  p-orbitals.  Isolated  S  atoms  also  are  predicted  to  bind  strongly  to  the  LB 
sites  of  Fe(110),  with  bridge  sites  found  to  be  transition  states  for  S  hopping  between 
neighboring  LB  sites.  MS  is  more  strongly  absorbed  than  FI2S  due  to  its  ability  to  form  a 
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primarily  ionic  bond,  and  S  is  the  most  strongly  bound  of  all  due  to  its  ability  to  form  two 
polar  covalent  bonds.  The  minimum  energy  paths  for  H2S  and  HS  dehydrogenation 
involve  rotating  an  H  atom  towards  a  nearby  surface  Fe  atom,  with  the  S-H  bonds 
breaking  on  the  top  of  only  one  Fe  atom.  The  barriers  to  break  the  first  and  second  S-H 
bonds  in  H2S  are  low:  0.1  eV  and  zero  on  Fe(110),  respectively,  suggesting  deposition 
of  S  on  Fe(110)  via  H2S  is  kinetically  and  thermodynamically  facile,  consistent  with 
experiments  [2]. 

We  then  considered  whether  alloying  the  steel  surface  might  be  able  to  inhibit  damaging 
surface  reactions  by  raising  the  barriers  to  molecular  dissociation  of  CO  and  H2S  [3]. 
We  chose  to  examine  Fe-AI  and  Fe-Si  alloys  because  they  are  the  only  binary  Fe  alloys 
with  favorable  thermo-mechanical  properties  that  might  allow  them  to  survive  under  the 
harsh  conditions  of  a  gun  barrel  firing.  We  had  already  determined  that  H2S  dissociates 
on  pure  iron  surfaces  much  more  easily  than  CO  does.  Although  FeAl  surfaces  raise  the 
barriers  for  H2S  dissociation,  they  significantly  lower  the  barriers  for  CO  dissociation.  On 
the  other  hand,  FesSi  surfaces  raise  the  barriers  for  CO  dissociation,  but  they  are  as 
vulnerable  as  Fe  surfaces  to  H2S  dissociation.  Our  findings  suggest  that  alloying  steel 
with  Al  or  Si  is  unlikely  to  simultaneously  increase  its  resistance  to  the  initial  stages  of 
chemical  degradation  by  CO  and  H2S.  However,  it  is  possible  that  alloying  the  outer  few 
microns  of  steel  with  Al  or  Si  still  might  prove  useful,  if  the  barriers  to  diffusion  of  C  and 
H  might  be  significantly  raised  by  alloying.  We  are  currently  investigating  this  possibility. 

The  fact  that  alloying  Fe  with  Al  lowered  dramatically  the  barrier  for  CO  dissociation 
made  us  wonder  about  the  possibility  of  an  industrial  spin-off,  namely  the  use  of  some 
FeAl  bimetallic  surface  as  a  potential  low  temperature  catalyst  for  CO  conversion.  We 
therefore  investigated  [4]  CO  adsorption,  diffusion,  and  dissociation  energetics  on  a 
monolayer  of  Al  covering  Fe(IOO)  [Al/Fe(100)].  We  predict  a  weakly  chemisorbed  state 
of  CO  to  exist  on  Al/Fe(100),  with  CO  adsorbing  on  the  fourfold  hollow  site  in  a  very 
tilted  fashion.  This  state  is  predicted  to  have  an  extremely  low  CO  stretching  frequency 
of  only  883  cm'\  indicating  a  dramatically  weakened  CO  bond  relative  to  gaseous  CO, 
even  though  the  molecule  is  predicted  to  bind  to  Al/Fe(100)  quite  weakly.  We  predict 
that  dissociation  of  CO  starting  from  this  weakly  adsorbed  state  has  a  barrier  of  only 
-0.35  eV,  which  is  -0.70  eV  lower  than  that  on  Fe(IOO).  To  understand  how  the 
underlying  substrate  changes  the  electronic  properties  of  the  supported  Al  monolayer, 
we  compare  CO  adsorption  on  Al/Fe(100)  to  its  adsorption  on  analogous  pure  AI(IOO) 
surfaces.  This  highly  activated  yet  weakly  bound  state  of  CO  suggests  that  Al/Fe(100) 
could  be  an  effective  low  temperature  bimetallic  catalyst  for  Fisher-Tropsch  type 
chemistry. 


B.  First  Principles  Exploration  of  Alternative  Protective  Coatings  for  Steel 

In  addition  to  examining  failure  mechanisms  for  gun  tube  erosion,  we  investigated 
possible  alternative  protective  coatings  for  gun  tubes.  Currently,  steel  gun  tubes  are 
coating  with  chromium,  which  has  two  problems.  One  is  that  the  Cr  layer  cracks  and 
reactive  gases  reach  the  steel,  leading  to  steel  erosion.  The  second  is  that  the 
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electroplating  process  that  produces  Cr  coating  on  steel  unfortunately  also  produces 
carcinogenic  Cr(VI),  an  environmental  hazard.  As  a  result,  the  Army  is  looking  for 
alternative  “green”  (environmentally  friendly)  coatings. 

Our  focus  for  coatings  has  been  on  ceramic  materials,  since  they  exhibit  low  thermal 
conductivity  and  high  melting  points,  both  of  which  are  crucial  in  the  gun  firing 
environment.  We  recently  examined  with  DFT  the  Cr/Fe  (100)  and  (110)  interfaces,  and 
predicted  that  the  ideal  work  of  adhesion  of  Cr  to  Fe  is  ~5.4  J/m^  for  both  interfaces, 
which  is  much  stronger  than  any  other  coating  we  have  examined  to  date  [5].  This 
anomalously  strong  adhesion  is  partly  due  to  the  lack  of  strain  (<1%),  but  is  primarily 
due  to  significant  spin  correlations  between  antiferromagnetic  Cr  and  ferromagnetic  Fe 
and  the  attendant  strong  d-d  bonding  that  results.  The  advantages  of  Cr  are  clear,  and 
will  be  hard  to  match. 

We  evaluated  MoSi2  as  a  possible  coating,  because  in  addition  to  favorable  thermal 
properties,  it  has  the  potential  to  be  corrosion-resistant  by  forming  a  silica  scale.  We 
found  that  MoSi2  strongly  adheres  to  Fe,  with  an  ideal  work  of  adhesion  of  ~3.9  J/m^  for 
two  low-index,  low-strain  interfaces.  This  value  will  be  a  lower  bound  to  measured 
adhesion  energies,  since  the  latter  will  be  larger  due  to  plasticity.  This  ideal  adhesion 
energy  for  a  ceramic  coating  to  Fe  is  much  stronger  than  predicted  earlier  by  us  for  ZrC 
and  TiC.  We  attribute  this  stronger  adhesion  to  increased  covalent  interfacial  bonding 
for  MoSi2/Fe  compared  to  metal  carbide/Fe  interfaces  (where  metallic  bonding  plays  a 
larger  role),  as  evidenced  by  the  rearrangement  of  electron  density  and  the  character  of 
the  local  density  of  states  upon  formation  of  the  interface  [6]. 

The  silica  scale  that  forms  on  MoSi2  should  help  protect  MoSi2  against  high-temperature 
corrosive  gases.  We  therefore  examined  the  interface  between  Si02  (in  the  cristabolite 
form,  used  as  a  model  for  amorphous  silica)  and  MoSi2  [7].  We  find  that  the  interfacial 
bonding  is  localized,  as  evidenced  by  an  ideal  adhesion  energy  that  changes  only 
slightly  with  the  thickness  of  the  Si02  layer.  Moreover,  the  adhesion  energy  displays  a 
relatively  large  (0.40  J/m^)  variation  with  the  relative  lateral  position  of  the  Si02  and 
MoSi2  lattices  due  to  changes  in  Si-0  bonding  across  the  interface.  The  most  stable 
interfacial  structure  yields  an  ideal  work  of  adhesion  of  5.0  J/m^,  indicating  extremely 
strong  adhesion.  Local  densities  of  states  and  electron  density  difference  plots 
demonstrate  that  the  interfacial  Si-0  bonds  are  covalent  in  character.  Mo-0  interactions 
are  not  found  at  the  Si02/MoSi2  interface  investigated  here.  Our  work  predicts  that  the 
Si02  scale  strongly  adheres  to  MoSi2,  and  further  supports  the  potential  of  MoSi2  as  a 
high-temperature  structural  material  and  coating.  Indeed,  it  is  important  to  create  this 
silica  scale  by  deliberate  oxidative  pretreatment  at  high  temperatures  (e.g.,  1200  °C),  in 
order  to  suppress  oxidatively-induced  fracture  of  MoSi2  that  may  occur  below  600  °C. 
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Development  of  the  following  algorithms  and  codes: 

•  Ultrasoft  spin-dependent  pseudopotentials  generated  on-the-fly  within  the  CASTEP 
planewave  DFT  code. 
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